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INTRODUCTION 


At Willow Run Laboratories in 1970, under NASA Contract 9-9784, a two 
channel analog ratio technique in the thermal IR wavelength region for map- 
ping gross compositional variations in silicate rocks was developed which 
Successfully produced automatic ratio recognition maps of quartz sand and 
sandstone near the Mill Creek, Oklahoma sand quarry. In the report which 
describes that effort [1] , a three channel IR technique was hypothesized 
for a more refined version of this ratio method, whereby the third channel 
is used for correction of temperature variations across the target scene. 

It was also hypothesized that the resulting temperature-corrected ratios 
would be at least crudely correlated with SK^ content of the geological 
targets encountered. 

During 1971, the state-of-the-art of geological remote sensing was 
further advanced by the processing of data from an aircraft flight. over 
Pisgah Crater, California. There are four major results from the Pisgah 
Crater data. First, the two-channel thermal IR analog ratio method from 
last year has been successfully used to map relative differences among var- 
ious silicate rock types. Secondly, the three-channel thermal IR tech- 
nique, hypothesized last year, has been implemented successfully to make 
digital temperature-corrected ratio maps which are more sensitive than the 
two-channel ratio maps to true compositional variations in geological tar- 
gets. Thirdly, a new technique involving the ratio of radiances in visible 
green to reflective IR channels has been theorized and qualitatively tested 
for the purpose of mapping variations in iron oxides on the surfaces of 
exposed rocks. 

The Pisgah Crater area in San Bernardino County was chosen as a test 
site for three reasons: there are many silicate rock types in close prox- 
imity, the semiarrid terrain is very sparsely vegetated, and considerable 
ground truth data was already available about the area from other reports 
and geologic maps [2,3]. The flight took place on 30 October 1971 between 
0800 (approximately 1 hour after sunrise) and 0840 hours, local time, at 
above-ground altitudes of 500 ft and 3000 ft. Only the 3000 ft data was 
processed through to the final stages because of monetary limitations. The 
morning was clear, with relative humidity less than 25% and a ground-level 
air temperature of approximately 5°C . This time of day was chosen because 
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temperature variations across the scene due to compositional effects (albe- 
do, thermal conductivity, etc.) were small, yet enough solar illumination 
was available for visible and reflective IR data collection. 

Dr. Larry C. Rowan, a U. S. Geological Survey geologist from the 
Regional Geophysics Office in Denver, and Ben Drake, a geology Ph.D. candi- 
date at The University of Michigan, greatly assisted the author in field 
checking these images. 

DESCRIPTION OF APPARATUS AND METHOD OF ANALYSIS 


The data were collected by The University of Michigan multispectral 
scanner (4) and a Honeywell two-element Hg:Cd:Te detector filtered (5) to 
the approximate 10% cutoff points of 8.2 pm - 10.9 pm and 9.4 pm - 12.1 pm, 
which will be called channels 1 and 2, respectively, aboard the University's 
C-47 aircraft. The two methods of analysis for the thermal IR data are 
described in references 1 and 5. Briefly, in the first method, the ratio 
of radiances in channel 1 to channel 2 was made the "signal" to the Uni- 
versity's analog computer. The resulting "ratio images" are bright where 
this R^ ratio is high and dark where R-^ is low. Silicate rocks display 
various spectral emittance minima in the 8 pm - 14 pm region, caused by 
reststrahlen bands (molecular vibrational modes) . Generally speaking, fel- 
sic (silica-rich) rocks tend to have spectral emittance minima at shorter 
wavelengths than mafic (silica-poor) rocks. Hence, there is a tendency 
for the R -^2 ratio to be lower for felsic rocks (dark in the ratio image) 
than for mafic rocks (bright in the ratio image) . 


The second thermal IR method involves the correction of R-^ for temp- 
erature variations across the scene. The spectral emittance of the geolog- 
ic target in a channel 3 (11.4 pm - 14.4 pm at 15% points) is assumed to 
be a constant value for all silicates, which is a reasonably good assump- 
tion in that wavelength region. From the channel 3 radiance, the rock 
temperature is determined (within 5°K) , and a ratio correction factor, Rg, 
is calculated. The R-^ ratio is then corrected for temperature variations 
across the scene to yield the temperature-corrected ratio 
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This three-channel scheme is done digitally. Laboratory emittance spectra 
and a theoretical atmospheric model are then used to calculate expected R 
ratios for various rock types. The measured values of R from the scanner 
data can be compared with these theoretical values to determine broad clas- 
sifications of rock type. A simpler two-channel scheme for temperature 
correction, which has proved qualitatively successful^ will be implemented 
quantitatively during the coming year. 


The method for the visible-reflective IR ratio images involves a sim- 
ple ratio of radiances in channel 5 (0.50 pm - 0.52 pm) to channel 7 (0.74 
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ym - 0.85 ym) ■ Electronic transitions of iron oxides produce a relatively 
large rise in spectral reflectance from 0.5 pm to .80 pm, compared with 
lesser or no rise in spectral reflectance over the same spectral region for 
felsic rocks ; and with a drop in reflectance for unoxidized mafic rocks. 
Thus, generally speaking, iron oxides appear dark, felsic rocks appear a 
medium gray tone, and unoxidized mafic rocks appear bright on the R, ratio 
image . 


RESULTS AND DISCUSSION 


A demonstration of the two-channel IR ratio method is given in figure 
1, which shows from top to bottom the analog images of channel 1 (8.2 ym - 
10.9 ym) , channel 2 (9.4 ym - 12.1 ym) , and R^ (the ratio of radiances in 
the two channels) for section 2A of the Pisgah data. This is a region 
approximately 5 miles south-southeast of Pisgah Crater. North is to the 
top of the images. The ratio image shows emittance variations which are 
indicative of rock type differences. The felsic mountains (D) and allu- 
vium (A) appear darker than the mafic lava (LA) , and the playa material is 
contrasted sharply against the alluvium. The topographically higher daci- 
tic mountains (D) are warmer than other parts of the scene, which makes 
the ratio for that region higher than it should be on the sole basis of 
reststrahlen position. The patchy appearance of the playa is primarily 
caused by emittance variations, because the gray level of the single chan- 
nel images are almost uniform across the floor of Lavic Lake, of which 
this is a part. The most important fact about this figure, however, is 
that felsic and mafic lithological units are clearly discriminated in the 
ratio image. 

The three-channel IR technique corrects the R^ ratio for temperature 
variations across the scene via equation (1) . It is demonstrated by the 
digital ratio map in figure 2, which includes almost all of the region 
shown in figure 1. The ratio R from equation (1) is divided into the fol- 
lowing 5 levels : 


Level 1 R < .969 

Level 2 .969 < R < .973 

Level 3 .973 < R 5 .981 

Level 4 .981 < R < .992 

Level 5 .992 < R 

The lowest recorded R in level 1 is .941 and the highest level recorded 
level 5 is 1.024. These values for R can be compared directly with the 
graph in figure 3> which plots R versus %Si02 as calculated from laboratory 
rock emittance spectra of Lyon [6] and an atmospheric model of Anding, et. 
al . [7], Levels 1 and 2 are primar-’ly felsic (acidic) silicates and levels 
4 and 5 are primarily mafic (basic) silicates. A graybody with constant 
emittance e = 0.94 throughout the 7 ym - 14 ym region would yield a ratio 
value of R = 0.975, which would appear in level 3, along with some of the 
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intermediate rocks. Graybodies of that type can be distinguished from sili- 
cates by ratioing the sum of radiances in channels 1 and 2 to the radiance 
in channel 3 for the purpose of determining the presence or absence of rest- 
strahlen bands, which was not done here. More and different ranges in R 
could have been printed, but these ranges and symbols were chosen for print- 
ing because they produced a more pleasing figure for demonstration purposes. 
The estimated error in R is ± .007. 

The primary effect of the temperature-correction can be seen in a com- 
parison of R for dacitic mountains (D) with R for the alluvium (A) in fig- 
ures 1 and 2. Figure 1, which is not temperature-corrected, shows the 
alluvium as having a lower R-^ ratio (darker in the R^ ratio image) than 
the adjacent dacitic mountains, even though the alluvium has a chemical 
composition similar to the mountains (the former consists primarily of rock 
fragments from the latter). However, the mountains, which are topographi- 
cally higher, were warmer than the alluvium at this time of morning (see 
the single channel images of figure 1), which biased the R-^ ratio toward 
a higher value in the mountains. Figure 2, on the other hand, shows nearly 
equal values for R,- the temperature-corrected ratio, in the alluvium and 
in the mountains. This is a good indication that the temperature- correc- 
tions have been successful. 

As a more stringent test of the quantitative accuracy of this method, 
the ratios of the dacite and plagioclase basalt samples of Lyon, plotted 
as data points in figure 3, can be compared with ratios for the dacitic 
mountains (D) (also for the dacitic fragments in the alluvium) and basalts 
of the Sunshine lava flows (LA) in figure 2. Lyon's dacite yields a theo- 
retical ratio of R = .947, which falls into level 1, a ratio frequently 
measured in the mountain (D) and alluvium (A) regions of figure 2 . Lyon's 
plagioclase basalt yields a ratio of R = 1.015, which falls into level 5, 
a commonly measured level in the Sunshine basalt flows (LA) of figure 2 . 
Therefore, the laboratory and scanner data quantitatively agree within 
reasonable limits (perhaps as small as the estimated ± .007 error in R) . . 

To demonstrate the third ratio technique, ratio images of radiances 
in the visible green, channel 5 (0.50 ym - 0.52 ym),to reflective IR chan- 
nel 7 (0.74 ym - 0.85 ym) were processed in a manner similar to the two- 
chanriel thermal ratio method discussed above) for approximately the same 
areas as shown in figures 1 and 2. Figure 4 shows, from top to bottom, 
analog scanner images of channel 5, channel 7, and the ratio R,--,. The 
ratio image shows gray level variations unlike the single channel images. 
Theoretically, iron oxide will appear dark in the ratio image. The two 
eruptive phases (labeled 1 and 2) of Sunshine lava display various shades 
of gray, which are darker than the alluvium (A) on the'west side. The 
younger phase 2 Sunshine lava appears darker in the ratio image than the 
older phase 1 lava. The greater presence of calcium carbonate on joint 
planes and in local minor depressions on the phase 1 Sunshine lava ac- 
counts, in part, for the lighter color of this unit in the single channel 
images. This may also partly account for the slightly brighter gray level . 



33-5 


in the ratio image. 

The alluvial deposits (A) on either side of the Sunshine lava appear 
quite different in the ratio image. It is probable that the drainage pat- 
terns carry more ferric oxides from the Pisgah and Sunshine lava flows to 
the eastern alluvium (surrounded by basalt) than to the western alluvium 
(basalt only on one side) . This would explain the darker gray level of the 
eastern deposits in the ratio image. However, the playa and alluvium east 
of the Sunshine flows appear to have similar spectral properties in these 
two channels of information, as evidenced by similar gray levels in the 
ratio image, which may indicate almost equal amounts of surficial ferric 
oxide in or on both types of material. 

Probably the most significant features in the ratio image of figure 4 , 
however, are the dark linear regions in the dacitic mountains on the left 
side of the image. The two prominent, dark, linear features running north- 
west to southeast are andesitic dikes in the predominantly dacitic mountain 
range. The longer dark feature, north of the alluvium inlet, is shown as 
an andesite dike on T. W, Dibblee's geologic map of the area [8]. The 
shorter dark feature is an andesitic dike which cuts across the dacitic 
peninsula on the southern border of the alluvium inlet. This short dike, 
which is not on Dibblee's map, may be a slightly offset continuation of the 
more prominent, longer dike. If so, there may be a west-to-east fault be- 
low the alluvium in the inlet. The andesitic dikes appear rusty brown to, 
the naked eye of a field geologist on the ground, and, if not in shadow, 
are easily recognized by ground observers. However, the dikes do not show 
up well in the single channel images. Contrarily, the ratio image produces 
excellent contrast between oxidized surfaces of andesite and the surround- 
ing dacite, whether in shadow or bright solar illumination. Finally, the' 
small dark spot at the top of the ratio image on the mountain-alluvium bor- 
der is a darkly weathered biotite quartz monzonite outcropping that shows 
up equally well in the ratio or single-channel images. 

Generally speaking, acidic (felsic) rocks should contain less ferric 
oxide than basic (mafic) rocks, because the latter have more iron available 
for oxidation. Thus, there should be some correlation between ar *d R57 

with respect to silicate rock type. This implies that, low R-^ and high 
R57 indicate acidic silicates, whereas high R]_2 and low R57 indicate basic 
silicates. Although desert varnish (ferric oxides) can cover any type of 
rock, it does not seem to prevent the acidic rocks from appearing brighter 
than basic rocks in the R^., ratio of figure 4 , except for cases where the 
felsic rocks contain appreciable amounts of biotite (biotite quartz monzo- 
nite in figure 4 ) . The weathered andesite in the dikes of figure 4 is an 
intermediate rock with relatively high iron content; hence, it has a high- 
er Fe2C>3 content and appears dark in the R57 ratio image, even though it is 
not bright in the Rq2 image of figure 1. 

There are several indications in the Pisgah data that the thermal in- 
frared R^2 ratio is not very affected by the amount of ferric oxides on 
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the rock surfaces. For instance, R 57 variations within Pisgah lava flows 
(not shown) are not accompanied by similar variations in Ri2- Laboratory 
tests are needed to confirm or deny this hypothesis. However, the iron 
oxide reststrahlen bands occur near 20 ym and greater, which hints that 
they may be relatively transparent in the silicate reststrahlen region be- 
tween 8 ym and 12 ym. 


CONCLUSIONS 

During the past year a 2-channel IR technique for discrimination among 
silicate rocks was tested for a second time, a three-channel IR method was 
tested for the first time, and a new visible-reflective IR ratio method was 
hypothesized and qualitatively tested for iron-oxide recognition. Both the 
two-channel and three-channel IR ratios methods were capable of discriminat- 
ing felsic from mafic rock types near Pisgah Crater, California. The three- 
channel ratio values measured by the scanner were found to agree, within 
reasonable limits, with the values of R calculated from laboratory data. 

The capability of the three-channel method to obtain absolute ratios is 
accompanied by a costly processing routine. The two-channel technique is 
superior for low thermal contrast scenes (AT <_ 5°C) about which some ground 
truth is available because of its speed and economy ($3400 for data process- 
ing of the first 150 data-miles and $750 for each additional 150 data-mile 
set; this does not include data collection costs). 

The reflective channel ratio method was found capable of enhancing the 
contrast of iron oxides in the presence of unoxidized or iron-deficient 
rock surfaces. Photographic IR methods can duplicate some of these visible- 
reflective IR ratio results [9, '10] with less image distortion and lower 
cost. However, this scanner ratio technique (which costs the same as the 
2 -channel thermal ratio processing) is superior for suppressing shadowing 
and atmospheric effects. 

All three techniques can be useful for geologic mapping. For example, 
the thermal ratio method should assist in the mapping of volcanic ash flows, 
which are sometimes indistinguishable to the naked eye from sedimentary 
material, and glacial tills, which can vary significantly in SiC >2 content * 
Similarly, the visible-reflective IR ratio should be useful in the explora- 
tion for gossans or othe r iron and nickel ore bodies, as well as for the 
mapping of iron-rich (lateritic) soils. These and other applications can 
be tested with esiting equipment and computer software. 
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FIGURE CAPTIONS 


FIGURE 1: 


FIGURE 2: 


FIGURE 3: 


Analog Infrared images of flightline 2 section A of a 
west-east flight about 5 miles south of Pisgah Crater. 
North is toward the top. Channel 1, channel 2, and the 
Rl 2 ratio image are at top, middle, and bottom, 
respectively. From left to right are dacitic mountains 
(D) , alluvium (A) , Sunshine lava flows (LA) , a 1 luvium 
(A), playa material (P) of Lavic Lake, and Pisgah lava 
(LA) of phase 1. 

Temperature-corrected digital recognition map of silicate 
rocks for most of flight line 2, section A (including 
most of the area covered in figure 1 ) about 5 miles south 
of Pisgah Crater. North is to the top. The top and 
bottom maps are contiguous. Felsic rocks (high XSiC^) 
appear yellow and mafic rocks (low %Si 02 ) appear blue. 

Temperature-corrected ratio R versus %Si02 for 25 rock 
samples of Lyon [ 6 ]. 


FIGURE 4: Analog visible and reflective IR images of flight line 2 

section A of a west-east flight about 5 miles south of 
Pisgah Crater (same area as in figures 1 and. 2). North 
is toward the top. Channel 5, channel 7, and R 57 ratio 
are at top, middle and bottom, respectively. From left 
to right are dacitic mountains (D) with dark (low ratio) 
andesitic dikes, alluvium (A), Sunshine lava (LA) of two 
eruptive phases, alluvium (A), playa material (P) , and 
phase 1 Pisgah lava (1) . 
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Temperature-Corrected Ratio RecognitionMap i 
of Silicate Rocks Near Pisgah Crater. Calif. Top 
and Bottom Maps are Contiguous Parts of Flight 
Line 2. Section A. (West is to the L6ft.) 


FIGURE 2 
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